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Is the silicate emission feature only influenced by grain size? 
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Abstract. The flattening of the 10 yum silicate emission feature observed in the spectra of T Tauri and Herbig Ae/Be stars 
is usually interpreted as an indicator of grain growth. We show in this paper that a similar behaviour of the feature shape 
occurs when the porosity of composite grains varies. We modelled the fluffy aggregates with inclusions of different sizes with 
multi-layered spheres consisting of amorphous carbon, amorphous silicate, and vacuum. We also found that the inclusion of 
crystalline silicates in composite porous particles can lead to a shift of the known resonances and production of new ones. 



1. Introduction 

The shape and strength of the silicate emission feature ob- 
served near 10 pm in the spectra of T Tauri and Herbig Ae/Be 
(HAeBe) stars is commonly used as a measure of grain growth 
in protoplanetary discs (see Natta et al. 2007] for a review). It is 
well-known theoretically that with an increase of the grain size, 
the feature becomes wider and eventually fades away. In the 
case of compact spherical grains with a composition of astro- 
nomical silicate, the 10 fim feature disappears when the grain 
radius exceeds ~ 2 //m (see Fig. |T}. 

The standard approach to modelling the 10 fim feature in- 
cludes calculations of light absorption by several populations 
of compact (and hollow) silicate spheres. We used amorphous 
and crystalline particles of small and large sizes to fit the ob- 
served emission profiles. The model was first suggested by 
Bouwman et al. (2001) and then modified by van Boekel et 
al. d20051 >. It was used by Schegerer et al. (120061 1 : Honda et 
al. (f2006l >: Kessler-Silacci et al. d2006b : Sargent et al. (f2006l >: 
Sicilia-Aguilar et al. (120071 1; and Bouwman et al. (120081 1 to fit 
the observational data. Further the authors searched for cor- 
relations between the estimated mass fractions of large and 
crystalline grains and different stellar and disc parameters like 
mass, luminosity, age, spectral type, etc. As a rule, the correla- 
tions are absent or rather weak (see Table 5, Sicilia-Aguilar et 
al.. 120071 

In this Letter, we show that the shape, position, and strength 
of the 10 fim silicate feature is also influenced by variations of 
the properties of small mass composite aggregate grains. We 
modelled the fluffy aggregates by multi-layered spheres (see 
also Voshchinnikov et al.. 120061 1 . This particle model permits us 
to include an arbitrary fraction of materials, and computations 
do not require large resources. 

2. Model of porous composite grains 



Fluffy particles should appear as a result of grain coagulation 
in interstellar clouds and protoplanetary discs (e.g., Henning 
& Stognienko, [T9961 Dominik & Tielens, [T9971 Jones, 120041 
Ormel et al., 2007). It is expected that aggregates can consist 
of several silicate and carbon sub-particles of different sizes 
that can be treated as inclusions in a vacuum matrix. 

We use the model of spherical multi-layered parti- 
cles introduced by Voshchinnikov & Mathis (1999). Later, 
Voshchinnikov et al. ((2005) demonstrated that the optical prop- 
erties of layered spheres resemble those of fluffy aggregates 
with inclusions of different sizes. 

Our model parameters are: the refractive indices and vol- 
ume fractions V,7V t otai of the materials and the radius of com- 
pact particles r compact . The amount of vacuum in the particle 
(the particle porosity P, < P < 1) is 



P = V vac /V to ,al = 1 - Vsolid/Vtotal, 



(1) 



where V so iid is the sum of the volumes of all species exclud- 
ing vacuum. The radius of porous particles is related to that of 
compact particles as 



'"compact '"compact 

rpc, u S = {l _ p)m = ( y solld /y toml) >/3- 



(2) 



The model of layered spheres combines all components, in- 
cluding vacuum in one particle (internal mixing) in contrast to 
the standard approach discussed above where "external mix- 
ing" (mixture of several individual populations of compact 
grains) is used. 

For calculations, we use different kinds of carbon and sil- 
icates: amorphous carbon Bel and AC1 (Rouleau & Martin, 
1991), amorphous silicate with olivine composition (Dorschner 
et al. , |19951 l, crystalline olivine (Fabian et al. . 120011 1. and astro- 
nomical silicate (astrosil; Laor & Draine, 1993). 
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3. Analysis of silicate emission 

The silicate feature in the N band was observed in spectra 
of a large variety of objects (see Henning 2008 for a recent 
review). We should caution the reader that we calculate ab- 
sorption efficiencies, but measure the fluxes. A detailed anal- 
ysis of disc spectra certainly requires radiative transfer calcu- 
lations. Assuming that the silicate emission is optically thin, 
we can compare observed fluxes with theoretical profiles. The 
latter are proportional to the product of particle absorption 
cross section by the Planck function with the particle tempera- 
ture oc C a bs(/l) B A (T,i). Radiative transfer calculations show that 
grains of different temperatures contribute to the silicate feature 
(see Fig. 1 in Schegerer et al., 12006). However, the dominant 
contribution for the 10 /mi feature comes from particles with 
Td = 200 - 400K for which the Planck function is approxi- 
mately constant in N band. Therefore, we can adopt that the 
shape of the feature depends primarily on the emission proper- 
ties of grainfl 

The profile of the feature can be described by the nor- 
malised absorption efficiency factor 2abs(^)/Gabs(^o) — 1, 
where the flux at wavelength Aq characterises the continuum. 
As usual, the value of Aq — 8.2 /mi is chosen (e.g., Schegerer et 
al.. l2006l 

Another representation of the optical behaviour is provided 
by the mass absorption coefficient, which is the ratio of ab- 
sorption cross section C a bs(A) to particle mass. In the case of a 
sphere of porosity P, it can be written as 

Qbsffl 3 Q abs (A) 
KsbsW - — — — - -; , (i) 

PdV total 4p djSO i id r, 

compact \ 1 * ) 

where pd = pd,soiid(l _ f) 1S the mean particle density. The 
value of pa, solid is obtained by averaging the density of all 
species excluding vacuum. In our calculations, we assume that 
Pd.si = 3.3g/cm 3 andpd,c = 1.85 g/cm 3 for silicate and carbon, 
respectively. 

Figure Q] shows the wavelength dependence of the nor- 
malised absorption efficiency factors for compact silicate 
spheres of diverse sizes. With a growth of the particle size (and 
mass), the 10 pm silicate feature broadens, its height decreases 
and the position of maximum shifts to longer wavelengths. 

In a similar manner, the silicate feature changes when the 
particle porosity grows (Fig. [2j(. However, in this case the parti- 
cle size increases only moderately, while particle mass remains 
the same. With a growth of porosity, the peak strength de- 
creases for normalised absorption (Fig. [2] upper panel) whereas 
the mass absorption coefficient becomes larger (Fig. 12 lower 
panel). The value of A^bs almost doubles at the peak position 
when we replace the compact particle by porous particle. 

It is well known that the shape and strength of the silicate 
feature depend on the type of the amorphous silicate, particle 
size, and fractal dimension (see, e.g., van Boekel et al., 2005; 
Schegerer et al., 120061 Min et al., 120061 1. Using the model of 
composite grains we can also investigate how carbon embed- 
ded in particles affects the characteristics of the silicate feature. 

1 Note that our case differs from the case analysed by Li et al. ( 2004 1 
where the cometary particles of different composition are located at 
the same distance from the Sun. 




Fig. 1. Wavelength dependence of the normalised absorption 
efficiency factor [(2abs(/l)/2abs(8.2 /mi)-l] for compact spheres 
consisting of astrosil. The effect of variation of the silicate 
emission shape with the particle size is illustrated. The obser- 
vational profile of the T Tauri starTr37 13-157 (Sicilia-Aguilar 
et al. , 120071 1 is shown for comparison. Note that this profile is 
not continuum subtracted. 
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Fig. 2. Wavelength dependence of the normalised absorption 
efficiency factor (upper panel) and mass absorption coefficient 
(lower panel) for layered spheres with r SjComp act — 0.1 pm. 
The particles consist of amorphous carbon (Bel), amorphous 
silicate with olivine composition (MgFeSi04), and vacuum. 
The volume fraction of components V,/Vtotai i s indicated. The 
particles are of the same mass (VBei/^soUd = VBei/(VBei + 
ViUgFeSiO-i) = 0-2), but of different porosity. The effect of varia- 
tion of the silicate emission shape with the particle porosity is 
illustrated. 
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Fig. 3. The same as in Fig. [2] (upper panel) but now for parti- 
cles of porosity P = 0.9. The particles are of the same porosity 
but of slightly different mass. The value of VBei/V so iid varies 
from 0.1 to 0.3. The effect of variation of the silicate emission 
shape with the volume fraction of carbon is illustrated. The ob- 
servational profile of the T Tauri star Tr 37 13-157 is shown for 
comparison. 

This dependence is plotted in Fig. [3] With addition of carbon, 
the feature very rapidly transforms into a plateau and then dis- 
appears. 

Comparing the data presented in Figs. [T]and[3] with obser- 
vations, it is possible to estimate characteristics of grains fit- 
ting the observations: r SiCompact * 2.0 fim or V B ei/Vsoiid « 0.2 
and r s , polous as 0.215 /mi (last value is obtained from Eq. © for 
P = 0.9). Hence the same observational data can be explained 
with particles whose radii differ by a factor of ~ 9 and masses 
by a factor of ~ 8800(1). 

Using the observations from Schegerer et al. (2006); van 
Boekel et al. (120051 ); Sargent et al. (12006b ; and Sicilia-Aguilar 
et al. (2007), we fitted observational and theoretical profiles of 
the silicate feature. We choose the best model by minimising 
the x 1 criterion. In all cases, the observational continuum is fit- 
ted by a straight line in the interval 8-12 fim for the sake sim- 
plicity. We refer the reader to Juhasz et al. (2008) for a detailed 
discussion of the physical models of continuum. We made the 
calculations for layered spheres with r s compact = 0.1 fim con- 
sisting of Bel and MgFeSiO^ 

The estimated porosity and amount of carbon in grains pro- 
ducing the silicate emission feature in protoplanetary discs are 
collected in online Table Q] At this stage, the stars with very 
pronounced crystalline peaks were eliminated from considera- 
tion. Table [T]includes 47 stars (30 T Tau and 17 HAeBe stars). 
The obtained values of P, VBei/Vsoiid, me rat i° °f masses of 
carbon to silicate M Be i /AfM g FeSio 4 and stellar age (if known) are 
given. Note that the determination of the age is often quite un- 
certain for pre-mainsequence stars (see discussion in Blondel 
& Tjin A Djie, 120061 1 . Therefore, we do not discuss possible 
correlations. 

The grain porosity exceeds 0.5 and the average value of P is 
equal to (P) = 0.64 + 0.15. Such particles resemble aggregates 
obtained both experimentally (Kempf et al., 1999) and theoreti- 
cally (Shen et al., 2008 0- The amount of carbon in grains is not 

2 See also discussion in Li et al. (2003 1. 
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Fig. 4. The observational normalised profile of the 10 /mi fea- 
ture for the Herbig Ae/Be star AB Aur (van Boekel et al., 2005 ; 
continuum is subtracted.) The best fit profiles are shown for 
particles containing amorphous carbon in the form Bel (dashed 
line) and AC1 (circles). The values of particle porosity P and 
fractional amount of carbon VBei/Vsoiid are indicated. 

very large (average volume fraction ( Vsei I V sou d) = 0.24 + 0.10 
and average mass ratio (MB e i/MM g FeSio 4 ) = 0.19±0. 12). These 
values increase by a factor of 3 - 4, if we replace Bel by the 
less absorbing amorphous carbon AC1. This fact is illustrated 
in Fig.|H where the results are given for both particle materials. 
Note that the particles containing AC1 are less porous. 

The variation of grain porosity without significant change 
of grain mass may explain the behaviour of silicate emission. 
This explanation is an alternative to the commonly used idea 
of large grains in protoplanetary discs. We will be able to de- 
cide between the two hypotheses after conducting spectropo- 
larimetry in the 10 /mi feature because a noticeable albedo 
of large grains manifests itself in polarisation of the scattered 
light. In this case, we expect unusual behaviour of polarisa- 
tion parameters (especially positional angle) within the feature 
profile in comparison with calculated profiles for dichroic ex- 
tinction (see, e.g., Henning & Stognienko, 1993; Prokopjeva & 
Il'in,|20D7). 

4. Crystalline silicates in composite grains 

Another interesting problem is the degree of crystallinity of 
dust in protoplanetary discs, which is related to the processes of 
partial grain evaporation and annealing (e.g., Gail, [2004). Due 
to the conversion of amorphous silicates to crystalline minerals, 
the particles may consist of different types of silicates. In order 
to show the effect of amorphous silicate matrix and vacuum 
on resonances produced by crystalline silicates, we calculated 
the feature profiles for composite particles containing Mg-rich 
crystalline olivine Mg, .gFeo.iSi04 as a component. The results 
are plotted in Fig. [5] where the upper panel illustrates the influ- 
ence of particle porosity on position and strength of emission 
peaks. It is seen that variations of spectra are significant: the 
shape of the feature changes (cf. Fig. [2), some peaks totally 
disappear and new peaks arise. A very pronounced peak with a 
maximum near A = 11.37 fim is observed for very porous parti- 
cles whereas compact and medium-porous particles have reso- 
nances near A = 10.89 /mi. The inclusion of crystalline silicates 



4 



Voshchinnikov and Henning: Is the silicate emission feature only influenced by grain size? 



O 



cryst. 



Mgi.gFe .iSi0 4 +vacuum : 

X X X X X X X 



lct -0.1/xni 




MgFeSi0 4 + cryst. Mg t 9 Fe iSi0 4 +vacuum 

. 4 v [orBterite /v„ lid = 




o 

w 

■a 2000 
1000 



Bel + am. MgFeSi0 4 + cryst. Mg 19 Fe 01 Si0 4 +vacuum 



XX X 



v for8tcrltc /v„ lld =o.o& 

' \ V B .i/W=0.30 




8 9 10 11 12 13 

A, /im 

Fig. 5. Wavelength dependence of the mass absorption coeffi- 
cient for compact and layered spheres with r SiCompact = 0.1 /im. 
The particles consist of crystalline olivine and vacuum (up- 
per panel), amorphous silicate, crystalline olivine, and vac- 
uum (middle panel) and amorphous carbon, amorphous sili- 
cate, crystalline olivine, and vacuum (lower panel). The vol- 
ume fraction of crystalline olivine and Bel is indicated. The 
position of peaks is marked. The thick solid line at the lower 
panel shows the mass absorption coefficient for the mixture of 
compact grains (external mixing) with volume fractions of con- 
stituents corresponding to solid curve (internal mixing). The ef- 
fect of variation of crystalline resonances with particle porosity 
and composition is illustrated. 

in a composite particle containing amorphous silicate (middle 
panel) or amorphous silicate and carbon (lower panel) changes 
the picture. The resonance near A * 11.4 /mi is clearly seen at 
the long-wavelength wing of the feature. Its position slightly 
shifts if the porosity changes. Fabian et al. ( 120011) found that 
such a peak appeared in the spectra of strongly-elongated par- 
ticles. A double peaked structure around A - 10.25 /im arises 
as well. This structure was not noticed previously in spectra 



of crystalline olivines (H. Mutschke, priv. commun., 2007). 
Note that, as expected, the mixture of separate constituents 
(amorphous silicate, carbon and crystalline silicate, thick line 
in lower panel) do not lead to the shift of peaks. Further cal- 
culations with different materials for wider wavelength range 
and detailed comparison with Spitzer observations (see Watson 
et al., 12007b might help to resolve the problem of grain crys- 
tallisation in protoplanetary discs and to answer the question 
whether the crystals occur in "isolation" or as part of porous 
grains. 
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Table 1. Grain porosity and fractional amount of carbon in grains as derived from fitting the 10 fim silicate emission feature. 
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h 


0.704 


0.149 


0.098 


3.3 


4 


age from (5) 


HD 36112 


h 


0.431 


0.122 


0.078 


4.4 


4 


age from (5) 


HD 37357 


h 


0.480 


0.206 


0.146 


10.0 


4 




HD 37806 


h 


0.732 


0.405 


0.382 


0.8 


4 




AB Aur 


h 


0.714 


0.187 


0.129 


4.8 


4 


age from (7) 


HK Ori 


h 


0.551 


0.400 


0.374 


5.8 


4 


age from (5) 


UXOri 


h 


0.392 


0.160 


0.107 


2.8 


4 


age from (5) 


V380 Ori 


h 


0.694 


0.493 


0.545 


7.4 


4 


age from (5) 



'Meaning of symbols: t = T Tau star, h = Herbig Ae/Be star. 

3 References for observational data: (1) Sicilia-Aguilar et al. (20071; (2) Schegerer et al. (2006 ); (3) Sargent et al. (2006 1; 
(4) van Boekel et al. (2005 1: (5) Blondel & Tjin A Djie < f2006t ; (6) Gudel et al. ( T2007T ); (7) Bertout et al. ( 120071 . 



